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ABSTRACT 

In this paper the analysis of reactive electrically conducting fluid flowing 
unsteadily through a channel formed by two infinite parallel plates, under the 
influence of a transversely imposed magnetic field with variable thermal 
conductivity is presented both ( for the linear and nonlinear equations ) 
analytically and numerically. The solutions for unsteady state temperature and 
velocity fields are obtained using unconditionally and convergent semi-implicit 
finite difference schemes. Also, to confirm the accuracy of the numerical 
solution, steady state solutions for velocity field and temperature field are 
obtained using regular perturbation series method. Results of the computations 
for velocity, temperature, and skin friction as well as Nusselt number are 
graphically displayed for various controlling parameters embedded in the 
problem. 
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INTRODUCTION 

Thermal properties, particularly thermal conductivity and diffusivity are essential materials 
parameters of bedrock controlling the heat transfer and temperature increases in the vicinity of 
repository. To accurately predict the flow and heat transfer rates, it is necessary to take into account 
the variation of thermal conductivity and viscosity (Elbarbary and Elgazery, 2004). Most of the 
practical situations demand for physical properties with variable characteristics. The materials 
properties such as viscosity and thermal conductivity are prone to vary with a temperature gradient 
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(Arunachalam and Rajappa, 1978). There has been considerable published work dealing with 
steady flow with variable thermal conductivity. Arunachalam and Rajappa (1978) determined 
forced convection in liquid metal with variable thermal conductivity and capacity in potential flow 
and derived explicit closed form of analytical solution. Chiam (1996, 1998) reported the effect of 
variable thermal conductivity on heat transfer assuming it to vary linearly with the temperature. 
Also mention may be made some of the recent works with variable thermal conductivity of (Subhas 
et al. (2009a), Subhas et al. (2009b) and Elbarbary and Elgazery (2004)). 

The research on unsteady natural convection problems has been continuously attracting the 
attention of many researchers due to it is wider physical applications in nuclear reactors, solid 
matrix heat exchangers, thermal insulation, oil recovery, surface catalysis of chemical reactions, 
dispersion of chemical contaminants in various processes, storage of nuclear waste materials and 
many others. Details on the applications of unsteady natural convection flow can be found in (Jha 
et al. 2011a). Recently Jha et al (2011a, 2011b, 2012 and 2013) reported an excellent works on 
unsteady flow of reactive fluid in a vertical channel as well as vertical tube. In all, the above- 
mentioned authors neglected the effect of reactive hydromagnetic fluid. Hydromagnetic reactive 
flows found applications in thermo-hydrodynamic lubrication of rotating thrust bearing (Das 
(1998) and Salem et al. (1980)). In this direction, Makinde and Beg (2010) presented the inherent 
irreversibility and thermal stability in a reactive electrically conducting fluid flowing steadily 
through a channel with isothermal walls. In another article Makinde and Chinyoka (2011) 
investigated numerically unsteady hydromagnetic generalizes Couette flow of a reactive third- 
grade fluid with asymmetric convective cooling. Recently Uwanta and Hamza (2014a) analyzed 
unsteady flow of reactive viscous, heat generating/absorbing fluid with Soret and variable thermal 
conductivity in a vertical channel formed by two infinite parallel plates. Uwanta and Hamza 
(2014b) studied unsteady natural convection flow of reactive hydromagnetic fluid in a moving 
vertical channel. 


The present work is committed to examine unsteady as well as steady natural convection flow of 
reactive hydromagnetic fluid with temperature dependent thermal conductivity in a channel 
formed by two infinite vertical parallel plates. 


Governing Equations 

We consider unsteady natural convection flow of reactive fluid, incompressible and electrically 
conducting under the influence of a transversely magnetic field of strength B 0 , having temperature 
dependent thermal conductivity in vertical channel formed by two infinite vertical parallel plates 
separated by a distance H as shown in figure 1 . Initially, when/ ' < 0 , the wall at y ' = 0, y ' = H 

and the ambient fluid are at the same temperature T 0 . At timet ' > 0 , the temperature of the wall 
y ' = 0 is instantaneously raised to T 0) and that of y ' = H is lowered to T 0 and thereafter 
maintained as it is, whereT^ >T 0 . We choose a Cartesian coordinate system with the x ' -axis along 
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the plate in the vertically upward direction and the y ' -axis is taken normal to the plate. All the 
fluid properties except density in the buoyancy term are considered as constant. Neglecting the 
reactant consumption (Makinde and Beg, 2010), using the Boussinesq’s approximation, we write 
the governing equations for the present physical problem in dimensional form as 


du ' _ c 

dt ' 3 

3 T' _ : 
3 1 ' pt 


T'=T n 


u' = 0 


B n 


T'=T 


u' = 0 


A 


A 


E' = C 


Reactive? 


Hydromagneticb 


Fluid 


y 


V 


y'=H 


The initial and boundary conditions for the present problem are 
t <0: w ' = 0, T ' = T 0 for 0 <y'<H' 

t >0: u' = 0,T ' =T a , at y ' = 0 > 

u' = 0,T'=T 0 , at y' = H 


( 1 ) 

( 2 ) 

( 3 ) 


Where a is the conductivity of the fluid, v is the kinematic viscosity, B 0 is the electromagnetic 
induction, j5 is the coefficient of volume expansion, Q is the heat of reaction, A is the rate constant, 
E is the activation energy, R is the universal gas constant, C 0 is the initial concentration of the 
reactant species, g is the gravitational force, C p is the specific heat at constant pressure, p is the 
density of the fluid, T 'is the dimensionless temperature of the fluid, T m is the temperature of the 
fluid at the plate , T 0 is the temperature of the fluid far away from the plate and K f is the thermal 
conductivity which is assumed to vary as a linear function of temperature in the form (Subhas et 
al. 2009a and Subhas et al. 2009b), K f (T ) = K 0 jl + m * {T '—T 0 ) j , where K 0 is the thermal 

conductivity of the ambient fluid and m* is a constant depending on the nature of the fluid. In 
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general, m* >0 for fluids such as water and air, while m* < 0 for fluids such as lubrication oils 
(Elbarbary and Elgazery, 2004). 


In order to solve the equations (1) to (3), we employ the following non-dimensional quantities and 
parameters 




t V 
H 2 


u 

— , Pr: 


ypc, 


,Gr 


e = E(T-T 0 ) l = QC () AEH 


RT n 


RT n 


exp 


f -E ^ 
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Evv 0 

§= mRTl_ 


,M 


0 T = 


GBlH 2 


vp 

E (R a> -T 0 ) „ RT { 


(4) 


RT„ 


,£ = ■ 


Using (4), the equations (1) to (3) can take the following form: 
du d 2 u 

— = — -+GrQ-Mu 
dt dy 


de_ 

dt 


(1 + 80) 3 2 # 5_\de_ 

Pr dy 2 Pr dy 


A 

-\ exp 

Pr 


0 ^ 
V l 4- CO y 


(5) 

( 6 ) 


With the initial and boundary conditions in dimensionless form are 


n = 0, # = 0, 0<y<l 
t>0:n=0, 6 = d T at y=0 > 
n = 0, 0 = 0, as y=l 


(7) 


Analytical solutions 

The governing equations (6) displayed in the previous section are highly nonlinear and exhibit no 
analytical solutions. The analytical solutions have played an important role in validating and 
exploring computer routines of complicated problems. Therefore we reduce the governing 
equations of this problem into a form that can be solved analytically by setting du/dt =0 and 
dO/dt = 0 to obtain the steady state version of the problem. Since at steady state no parameter as 
any significant effect, the steady state solutions of (5) and (6) with the help of (7), using 

perturbation series method of this form u — ^ Uj A J and 6 - ^ 0 j A 1 for velocity and temperature 

7=0 y=0 

respectively by taking 0 T = 1 can be written as 
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u = 


A cos U^M y ) +B sinh( 4m y y 

M M 


+; i 

+A 2 


C cosh (Vm y)+D sinh (4m y ) +H l y 4 +H 2 y 3 +H 3 y 2 +H 4 y +H 5 
E x cosh (\[m y)+E 2 sinh (4m y ) J rL l y 1 +L 2 y 6 +L 3 y 5 +L 4 y 4 +L 5 y 3 + L 6 y 2 +L 7 y +L 8 


( 8 ) 
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(9) 


Using (8) and (9), we write the steady state skin friction and Nusselt number on the lower boundary 
as 
du 


dy 

90 

dy 


y=0 


f \ 

b 4m~ a(d 4W +h 4 )+a 2 (e 2 4W + l 7 ) 


— — 1 + Ak , + A 2 k * 
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Numerical Solutions 

The coupled nonlinear unsteady partial differential equations (5) and (6) are solved numerically 
using a semi-implicit finite difference scheme given by (Makinde and Chinyoka, 2011). We used 
forward difference formulas for all time derivatives and approximate both the second and first 
spatial derivatives with second order central differences. The semi-implicit finite difference 
equation corresponding to equations (5) and (6) are as follows 




+ (l + 2r l )u < j N+1) -r x u^ l) = r 2 u\ N _t + (l-2r 2 - AtM )u\ N) + r 2 u ( /[l + AtGrO 


-r^ + (Pr+ 2 r,)6f +1) - +1) 


'j + 1 


(r, + r, )«>",' + (Pr- 2r, - 2r, )«<" 1 + (r, + r, )«",> + 


(C-Cf+^ e xP 


Q(N ) \ 

\ + ee {N ) , 

V j J 


( 12 ) 

(13) 


where r, = gAt / Ay 2 , r 2 = (1 - g)At / Ay 2 , r 3 = SAt 6" /Ay 2 , r 4 = SAt / A Ay 2 and 0 < % < 1 . We 
chose <^ = 1 , the detailed reasons to this particular selection is documented in [14]. Also the 
analytical solutions displayed in the previous section are used as a check on the accuracy and 
effectiveness of the numerical scheme. Again, in order to reconfirm the accuracy of the scheme, 
the numerical results for velocity and temperature are compared with the analytical solutions. At 
larger time a steady state condition is reached as can be seen in figure (2), an excellent agreement 
is reached between these results. 



Fig. 2. Unsteady and steady state solutions for velocity and temperature profiles 
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RESULTS AND DISCUSSION 

The numerical investigations were carried out for positive values of the buoyancy parameter 
Gr > 0 (which corresponds to the cooling problems). The cooling problem is often encountered in 
engineering applications such as the cooling of nuclear reactors and cooling of electronic 
components. Unless otherwise stated, the values: X = 0.1, Gr = 1, M = 1 , 0 T =1, Pr = 0.71, 

c> = 0.1 and £ - 0.01 are used for the investigation, while others are varied in some instances. 
Results obtained are presented graphically for velocity, temperature, skin friction and Nusselt 
number for various flow parameters. 


Fig-3. Variation of unsteady and steady state temperature for X and t 


Fig-4. Variation of unsteady and steady state velocity for X and t 

Figures 3 and 4 shows the effects of Frank- Kamenetskii parameter and dimensionless time on both 
the temperature and velocity respectively .These figures disclose that as time increase, the 
temperature and velocity increases until a steady state is reached. It is also observed that as X 
increase, both the temperature and velocity profiles increases. These harmonize with the results 
reported by (Jha et al. 2012). 
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Fig.5. Variation of unsteady and steady state temperature for 5 and t 



Fig.6. Variation of unsteady and steady state velocity for S and t 


Figures 5 and 6 illustrated the temperature and velocity distributions of fluid for different values 
of the variable thermal conductivity and non-dimensional time. It is noted that both the temperature 
and velocity increases with increasing non-dimensional time until a steady state is reached. It is 
also observed that both temperature and velocity increases with increasing values of the variable 
thermal conductivity parameter. 



Fig.7. Variation of unsteady and steady state velocity for Gr and t 
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Fig.8. Variation of unsteady and steady state velocity for M and t 


Figures 7 and 8 show the velocity profiles of the fluid for different values of the buoyancy and 
magnetic parameters as well as the non-dimensional time respectively. In figure 7 it is found that 
the velocity increases with increasing dimensionless time but reverse effect is observed in figure 
8. Also in figure 7 it is seen that increasing value of Gr imply an increase in the buoyancy forces 
relative to the viscous forces and this reflected in the progressive increase in the vertical velocity 
of the flow. In figure 8 it is observed that increasing magnetic parameter decreases the velocity 
throughout of the flow region bounded by the channel walls. The physics behind this happening is 
that the presence of transversely applied magnetic field produces a hold back force (Lorentz force) 
similar to the drag force, which shrink the velocity. 



Fig.9. Variation of unsteady and steady state temperature for Pr and t 


Figure 9 displayed the effect of the Prandtl number on temperature and dimensionless time. The 
Prandtl number was chosen Pr = 0.71, 1 and 7 which corresponds to air, electrolyte solution and 
water respectively. It is observed that as the time increase and the Prandtl number increase the 
temperature of the fluid decreases. 
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Fig. 10. Variation of unsteady and steady state skin friction at y = 0 against A 



Fig.ll. Variation of unsteady and steady state Nusselt number at y = 0 against A 


Figures 10 and 11 illustrated the skin friction and Nusselt number for the different values of A 
and dimensionless time. In figure 10 and 11 it is seen that as the time increases, the skin friction 
and Nusselt number increases until a steady state condition is reached. Also from these figures it 
is observed that as A increases both skin friction and Nusselt number increases. 



Fig.12. Variation of unsteady and steady state skin friction at y =0 against S 
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Fig.13. Variation of unsteady and steady state Nusselt number at y - 0 against S 


Figures 12 and 13 displayed the skin friction and Nusselt number for different values of non- 
dimensional time and the variable thermal conductivity parameter. From these figures it is noted 
that both skin friction and Nusselt number increases by increasing the non-dimensional time until 
a steady state is attained. It is also observed from figure 12 that increases in the variable thermal 
conductivity parameter have no significant effect on the skin friction. From figure 1 3 it is seen that 
increases in the variable thermal conductivity leads to the increase in the Nusselt number. 



Fig. 14. Variation of unsteady and steady state skin friction at y - 0 against Gr 



M 


Fig. 15. Variation of unsteady and steady state skin friction at y - 0 against M 
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Figure 14 and 15 represent the skin friction for different values of buoyancy, magnetic parameter 
and dimensionless time. From these figure it is observed that the skin friction increases by 
increasing the non-dimensional time until a steady state is attained. In figure 14 it is also seen that 
the skin friction increases by increasing the buoyancy parameter. From figure 15 higher values of 
magnetic parameter tend to reduce the skin friction 


CONCLUSION 

This paper studied the problem of unsteady natural convection flow of a reactive hydromagnetic 
fluid with variable thermal conductivity. The temperature and velocity field for steady state version 
of the problem are obtained analytically using perturbation series method. The complete forms of 
the problem are solved numerically using unconditionally and convergent semi-implicit finite 
difference scheme. It is observed that both temperature and velocity distributions increases with 
increasing non-dimensional time until a steady state is reached. It is also further noticed that as the 
reaction strength and temperature dependent thermal conductivity increases, the velocity and 
temperature increases as well. Higher values of magnetic parameter retard the motion of the fluid, 
while greater values of buoyancy parameter accelerate the motion of the fluid. 
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